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Abstract: To predict the fractal decoded image quality more efficiently, an accumulated collage error 
coefficient (ACEC) based method was proposed in this study. Firstly, the definition of ACEC was 
introduced to describe the relationship among the upper bound, lower bound, and actual value of the 
accumulated collage error of all range blocks. Moreover, the definition and monotonicity of the relative 
error of ACEC were also defined and discussed. While the relative error of ACEC reaches a relatively 
small value, the average collage error (ACER) can be estimated approximately, and then the encoding 
process can be terminated to directly predict the peak signal-to-noise ratio (PSNR) quality of decoded 
images. Experimental results show that compared with the previous method, the proposed method can 
predict the decoded image quality more accurately with fewer computations. 
Keywords: Fractal image coding; Decoded image quality; Accumulated collage error coefficient 
1. Introduction 

Images play the most important role in human perception, and everyday a huge number of 
images are stored, processed, and transmitted, which will lead to great pressure to the 
transmission and storage requirements. Efficient image compression technology is an 
effective way to address the problem of reducing the amount of data. Different from 


conventional image compression techniques such as discrete cosine transform and discrete 


wavelet transform, fractal image coding realizes image compression by establishing an 
iterated function system (IFS) whose fixed point can well approximate the original image. 
The kernel idea of fractal image coding was firstly introduced by Barnsley, and then Jacquin 
proposed the first block based fractal coding method!'?!, Although fractal image coding has 
the advantages of novel idea, fast decoding, potential high compression ratio, and resolution 
independence!?“!, the drawback of fractal image coding was the high computational 
complexity in encoding process. Thus, fast fractal image encoding becomes one of the key 
research directions. In recent years, many fast fractal encoding methods were proposed by 
converting global block matching into local block matching'-!. For example, Chaurasia 
proposed composite statistical features, and Gupta extracted the features of the blocks in the 
DCT domain. The above two methods located the nearest neighbors in the feature space as 
the candidates of best-matched domain blocks, and then the block-matching operations were 
only carried out within the candidates. Instead of finding out the best matched domain block, 
Zheng specified the update times of fractal codes to obtain an acceptable domain block, and 
the encoding process can be effectively shortened. Further, the no-search fractal image coding 
methods directly assigned one domain block as the best-matched domain block for each 
range block without searching operations, and thus real-time encoding and higher 
compression ratio can be realized at the expense of poor decoded image quality!*!!!, After 
many years’ development, fractal image coding has been gradually applied to other image 
processing applications, such as image denoising!'*"'*! , image retrieval!'!??!!, image super 
resolution!*?6, watermarking?’**!, image hashing®°?"!, and head pose estimation??~91, 


The fractal decoded image quality prediction method was a recently proposed technique 


which can predict the decoded image quality only with partial encoding process4!. In the 
previous method, the average percentage of accumulated collage error is used to predict the 
decoded image quality, but the proportional relationship among the upper bound, lower 
bound, and actual value of accumulated collage error has not been considered. Thus, to 
predict the decoded image quality more efficiently, an accumulated collage error coefficient 
(ACEC) based method was proposed in this study. Based on the upper and lower bounds of 
the accumulated collage error of all range blocks, ACE%;, and ACE‘, , the definition of 
ACEC was proposed to describe the proportional relationship among ACE%,, ACEX,,, and 
the actual value of the accumulated collage error of all range blocks. Further, the definition 
and monotonicity of the relative error of ACEC were also defined and discussed. In encoding 
process, while the relative error of ACEC reaches a relatively small value, the encoding 
process can be terminated, and the average collage error (ACER) can be approximately 
estimated. Finally, based on the logarithmic relationship between ACER and the peak 
signal-to-noise ratio (PSNR) quality of decoded images, the decoded image quality can be 
directly predicted. Experimental results show that compared with the previous method, the 
proposed method can predict the PSNR quality of decoded images more accurately with 
fewer computations. The main contributions of this study can be summarized as: 

1. Based on the upper and lower bounds of the accumulated collage error, the definition of 
accumulated collage error coefficient (ACEC) was proposed to describe the relationship 
among the upper bound, lower bound, and actual value of the accumulated collage error of all 
range blocks. 


2. How to estimate ACEC, the definition and monotonicity of the relative error of ACEC 


were also introduced and discussed, based on which the proposed ACEC based prediction 
method can provide higher prediction accuracy with fewer computations than the previous 
method. 

3. The proposed method can be combined with the existing fractal encoding methods to 
predict the decoded image quality only with partial encoding process. 

This paper is organized as follows: Conventional fractal image coding is reviewed in 
Section 2. The principle and detailed procedures of the proposed method are described in 
Section 3. The experimental results are presented and analyzed in Sections 4 and 5, 
respectively. The conclusion is given in Section 6. 

2. Conventional fractal image coding 

Fractal image coding aims to establish an IFS whose fixed point can approximately 
approach the original image. The first block-based fractal encoding method was proposed by 
Jacquin, and the detailed procedures can be described as!!!; 

Step 1: Divide the M x N input image f uniformly into nonoverlapping B x B range blocks 
R;, i=1,2,3,....NumR, where NumR denotes the total number of range blocks. 

Step 2: Slide a 2B x 2B window over the original image horizontally and vertically, and the 
domain blocks Dj, j=1,2,3,....NumD, can be obtained, where NumD denotes the total number 
of domain blocks. 

Step 3: Establish a domain block pool (DBP) by contracting the domain blocks to the same 
size of range blocks. Then, an extended domain block pool (EDBP) can be obtained by 
performing eight isometric transformations on each domain block within DBP. 


Step 4: The best matched domain block for each range block can be obtained by searching 


within EDBP and minimizing the following function: 


CE(R,) = min (min |R -s,D,- oaf) 


(1) 
i=1,2,3,..., NumR, j =1,2,3,...,. NumDx8 


where CE(Ri) denotes the collage error of R;. I denotes a B x B matrix whose components are 
all ones. s; and o; denote the scaling and offset coefficients of the affine transformation, 
respectively. D; is the best matched domain block of Ri. Further, s; and o; can be computed 


with the least square method as: 


s, = (RFI, D,-d,1)/|p,-d,1| » 0, =F sd, 
i=1,2,3,..., NumR, j =1,2,3,..., NumDx8 


(2) 


where (**) denotes the inner product. 7 and d, denote the averages of R; and D,, 
respectively. 

In fractal decoding process, arbitrary M x N image can be selected as the initial image. 
Then, the same transformations as those established in encoding process are performed to 
reconstruct all range blocks until one iteration is completed. After about 10 iterations, the 
decoded image can be obtained. Figure 1 (a) illustrates the rotated 256 x 256 Bird image as 
the initial image, and Figure 2(b)-(f) illustrate the first five iteration images in decoding 


process while the Bird image is encoded in encoding process. 


(a) Initial image (b) First iteration (c) Second iteration 


(d) Third iteration (e) Fourth iteration (f) Fifth iteration 


Fig.1 Initial and the first five iteration images in decoding process. 


3. Prediction of the decoded image quality with the accumulated collage error 
coefficient. 
3.1 Definition of the accumulated collage error coefficient 


Substituting Eq.(2) back to (1) yields: 


CE(R,) =—|R -s,D,-o,1| = (ir -7I[ -s?|D, -arf | 
-FID -d IV (3) 
=L |R -7/ Ke ce aay = L JR -7I (I-LCc?) 
p |R -72| |D, -4,1| p 


where LCC, denote the linear correlation coefficient between R, and D, . Because 
0<LCC <1, we have 


1 = 
CE(R,) S elk -rI 
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(4) 


From Eq.(4), we know that for arbitrary range block, the variance provides the upper 
bound for the associated collage error, and only the range blocks with large variances may 
have large collage errors. Thus, all range blocks are sorted by their variances from largest to 
smallest and encoded block by block. In this process, the accumulated collage error (ACE) of 


all range blocks can be described as: 


EM! = -5 CE(R) ==> 


- 7a [1-(ecc,y | 6) 


Then, if we use LCCun and LCC, to represent the minimum and maximum of all 


LCCs, the upper and lower bounds of ACE*“" can be described as: 


ACE! = S5 CE(R = 
Max ( Max 
i=l 


raf [1-(L0C w) 


NumR NumR (6) 


ACEi = 2, CE(R,) sin a S |R, -7 Al [1-(LCCpax j] 


i=l 


In encoding process, the range blocks can be divided into two categories: Coded range 


blocks and uncoded range blocks. The ACE of coded range blocks can be represented as: 


Num 


ACE = È cee) = > ao ral |) nce) | (7) 


Coded 


where Num denote the number of coded range blocks. Similarly, the ACE of uncoded 


range blocks can be represented as: 


Num Uncoded 1% coded 


ACE = = y, CE(RO™*) = = =: |i = FI 2 [i E [ECE y (8) 
i=l E 


Uncoded Uncoded n d 


where Num denote the number of uncoded range blocks. If we use LCCy,,. 
LCCH! to represent the maximum and minimum of the LCCs of uncoded range blocks, 


the upper and lower bounds of the ACE of uncoded range blocks, ACEy<" and 


Uncoded 
ACE mi ~ , can be expressed as: 


ACEO = ei CER). = È |Re -7af i eee) | 
ACE = ee CERT) 2 3= 2 S S jeee- rl i (LOC f| 7 
Then, we have the inequality as: 
ACEi SAGE.” <ACER (10) 


Moreover, because the uncoded range blocks are the subset of all range blocks, LCCyf*""" 


and LCC" can be approximated by LCC,,,, and LCCy,;,, respectively. Then, Eq.(9) 
pp y 
can be rewritten by replacing LCCy<" and LCCY" with LCCy,, and LCCyuin , 


respectively, as: 


ACE Upota = Bs CERI") = + > Revere = ral [i = (eG A 

Uncoded umUtoded (11) 
ACE Vneoded 2E 5 ` CER?) = = > arr wes rl [i = (EOU )'| 

i=l i=l 


For all range blocks, to describe the proportional relationship among ACE“, ACE*.. 
and ACES , the ratio of ACE, to ACE" ratioto ACE*" to ACE), ratio is defined 


as the accumulated collage error coefficient, ACEC*", which can be calculated by: 


a _ ACES! /ACE“" ACES! ACES! 
ACEC™ = All Al 2 (12) 
ACE™/ACE vin (ACE4") 
Further, Eq.(12) can be expanded as: 
All All 
ACEC*"= ACE max ACE min (13) 


(ACES + ACE Unceded j 
Further, by Eq.(10), the upper and lower bounds of ACEC“!", ACEC“! and ACEC%,. , 


can be represented, respectively, as: 


ei, Rees 
Max (ACE AACE ord j 
Min 
(14) 
ACEC! = ACE vex ACE sin 


Min ACES 4. ACR Uncoded j 


Max 


While ACEC}, can approach ACEC}\, sufficiently, we can estimate ACEC" as: 


pa All All 
ACEC" f ACEC Max A (15) 


3.2 Relative error of ACEC and Predicting the PSNR quality of decoded images 


Based on the dynamic range of ACEC^" , i.e., the deviation between ACEC! and 


ACEC", the relative error of ACEC*", REaczc, can be defined as the ratio of the half 


deviation between ACEC, and ACEC}; to the estimated ACEC“" by Eq.(15) as: 


(ACEC -ACEC } /2 


Max 


ACEC — rT (16) 
ACEC 
By Eqs.(15), (16) and (14), REacec can be expanded as: 
(ACEC i, = ACEC (a, )/2 sACECM. — ACEC*! 
i ACEC" ACEC, ACEC, 
17 
(AcE +ACE ei j 2 (ACES + ACE noted j ( ) 


Max 


( ACEC" 4 ACE Yri j F ( ACEC! + ACE j 


Because a sufficiently large number of domain blocks within EDBP can guarantee that 


LCCua 71, ACERO"! ~ 0 , and Eq.(17) can be rewritten as: 


>’ Min 


(ACES + ACE Ye i _ (ACES j 


ACEC ~ ( ACE! 4. ACE Unodd ) re ( ACEC j (18) 


Max 


By setting the derivatives of RE crc regarding ACES and ACERS"! to zeros, 


respectively, we have 


Eo =-4 (ACE Coded ) (ACE eee ) (ACE Coded 4 ACE Cee )< 0 
(19) 
a | (ACES j (AcE + ACH Woes ) $0 
Max 


At the early stage of encoding process, ACES is much smaller than ACEy<°"" . By 
Eq.(18), we know that RE, =! . As the encoding process proceeds, because 


ORE ycrc aa ORE cec 


GACE C4 <0 CACE Ydd > > RE sane decreases monotonically while 


ACE“ increases, and ACE‘ decreases, respectively. At the final stage of the encoding 
process, ACES is much larger than ACERS"? . By Eq.(18), we know that RE crc 0. 


Thus, REacec set at a relatively small value can guarantee that a sufficiently large number 


of range blocks have been encoded, and thus the maximum and minimum of the LCCs of 
coded range blocks, LCC" and LCC? , can be used to approximate those of all range 
blocks, LCCua and LCC \,, respectively, as: 


LCCC ~ LCC 
| Max Max (20) 


LEC E LOC 
By Eq.(20), Eq.(6) can be rewritten by replacing LCCy,, and LCCy;, with LCC" 


and LCCH , respectively, as: 


ACE“! = YT CER) -1 S$" |r -Fll j1-(Lecgae) | 
Max = i 7 Max B? i i Min 


i=l 


NumR ] Numer b ; (21) 
ACES = 2, CER )w = z7 2, |R - Fl] 1-(Locg) | 
i=l i=l 
and similarly, Eq.(11) can be rewritten as: 

ACE = uu CE( R ass pe 1 MS Roe = i; I 2 j= Lec 2] 

i=l B i=1 

p an _ (22) 
ACE — > CE( RSS ee = ayn Roe was A I 2 1 a LCC. 2 

= Be £ J 


Thus, ACEC and ACEC in Eq.(14) can be obtained. Moreover, by Eq.(16), small 
RE csc can also guarantee that ACEC). is relatively close to ACEC). Then, ACEC™ 


can be estimated by Eq.(15), and ACE*" can be computed by Eqs.(12) and (21) as: 


nce"= [ACEACES 


AI (23) 
ACEC 
Further, the average collage error (ACER) of all range blocks can be computed as: 
ACE 
ACER = (24) 
NumR 


In addition, the peak signal-to-noise ratio (PSNR) was used to measure the decoded image 


quality as: 
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PSNR = log, fass/ a 5 >>I f(b j- PO (i, in| (25) 
jal il 
where f and f°" represent the input and decoded images, respectively. Moreover, 
there exists a logarithmic relationship between ACER and PSNR as: 

PSNR = £, + B,log,, (ACER ) (26) 
where /, and £, are both constant values. Thus, by Eqs.(23) and (24), we can directly 
predict the PSNR quality of the decoded image. The detailed procedures of the proposed 
method are listed as follows: 

Step 1: Divide the input image into range and domain blocks, and sort the range blocks 
by their variances from largest to smallest. 

Step 2: Take out one uncoded range block and encode it. Based on Loci" and 
Loci | ACECN, and ACEC}, can be calculated by Eqs.(14), (21) and (22). Further, 
RE xcec can be computed by Eq.(16). If RE,cuc is smaller than 0.15, go to Step 3. If not, 
turn back to Step 2. 

Step 3: Calculate ACER by Eqs.(24), (23), (21), (15), and (14), and then predict the 
decoded image quality by Eq.(26). 


4. Experiments 


(a)Bird (b) Owl (c) Trees 
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(d) City (e) Cat 


Fig.2 Five test images. 

Five 256 x 256 images (in Fig.2), Bird, Owl, Trees, City, and Cat, are adopted in our 
experiment. Three state of the art fractal coding methods and Jacquin’s method are used to 
assess the performance of the proposed method. The scaling and offset coefficients, s and o, 
are quantized by 5 and 7 bits, respectively. The previous prediction method is compared with 
the proposed method by prediction accuracy and the number of the computations required. 
Table 1, 2, 3, and 4 list the experimental results of Jacquin’s, Chaurasia’s, Zheng’s, and 
Gupta’s methods, respectively. In Tab.1, while the size of range blocks is 4 x 4, for Jacquin’s 
method, the PSNRs of decoded images are listed in the second row. Moreover, in the third 
row, 4096 range blocks are needed to be encoded for all test images, and the percentages of 
the computations required (PCRs) are all considered as 100%. For the previous method, the 
fourth to sixth rows list the predicted PSNRs, deviations between the predicted and actual 
PSNRs, and PCRs. Similarly, the seventh to ninth rows show the counterparts for the 
proposed method. From Tab.1, we see that if the size of range blocks is 4 x 4, the average 
deviation decreases from 0.03dB to 0.02dB, the average PCR decreases from 62.72% and 
56.20%. If the size of range blocks is 8 x 8, the average deviations decreases from 0.15dB to 


0.12dB, the average PCR decreases from 62.81% and 54.47%. Thus, we know that for 
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Jacquin’s method, the proposed method can provide higher prediction accuracy with fewer 
computations than the previous method. 

Further, Table 2, 3, and 4 list the experimental results of Chaurasia’s, Zheng’s, and Gupta’s 
methods, respectively. In Table 2, while the sizes of range blocks are set at 4 x 4 and 8 x 8, 
the average PCRs decrease from 62.98% to 56.46% and from 62.91% to 54.65%, respectively, 
and the average deviations decrease from 0.09dB to 0.07dB and from 0.18dB to 0.09dB, 
respectively. In Table 3, when the sizes of range blocks are set at 4 x 4 and 8 x 8, the average 
PCRs decrease from 70.97% to 64.64% and from 73.49% to 64.17%, respectively, and the 
average deviations remain at 0.05dB and decrease from 0.18dB to 0.10dB, respectively. In 
Table 4, when the sizes of range blocks are set at 4 x 4 and 8 x 8, the average PCRs decrease 
from 63.54% to 56.97% and from 63.11% to 54.86%, respectively, and the average deviations 
decrease from 0.08dB to 0.06dB and from 0.17dB to 0.11dB, respectively. In summary, the 
proposed method can predict the PSNR quality of decoded images more accurately with 
fewer computations than the previous method. 


Tab.1 Performance comparison between the previous and proposed methods for Jacquin’s method"), 


Test images Bird Owl Trees City Cat Average 
. Decoded(dB) 31.68 27.25 33.89 31.86 29.45 x 
Jacquin’s!!! ; 
Computations(%) 100 100 100 100 100 x 
4 Predicted(dB) 31.66 27.22 33.80 31.87 29.43 x 
Previous41 Deviation(dB) 0.02 0.03 0.09 0.01 0.02 0.03 
4 Computations(%) 68.99 72.73 71.17 50.90 49.80 62.72 
Predicted(dB) 31.67 27.24 33.82 31.87 29.45 x 
Proposed Deviation(dB) 0.01 0.01 0.07 0.01 0.00 0.02 Y 
Computations(%) 60.77 66.24 64.31 44.38 45.29 56.20 Y 
f Decoded(dB) 26.74 22.59 28.78 27.38 25.31 x 
Jacquin’s!!! 
8 Computations(%) 100 100 100 100 100 x 
x Predicted(dB) 26.63 22.37 28.62 27.30 25.11 x 
8 Previous’ Deviation(dB) 0.11 0.22 0.16 0.08 0.20 0.15 


Computations(%) 67.29 74.22 66.21 56.25 50.10 62.81 
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Predicted(dB) 26.65 22.40 28.69 27.31 25.16 x 
Proposed Deviation(dB) 0.09 0.19 0.09 0.07 0.15 0.12V 
Computations(%) 55.47 66.80 57.13 49.32 43.65 54.47 Y 


Tab.2 Performance comparison between the previous and proposed methods for Chaurasia’s 


method", 
Test images Bird Owl Trees City Cat Average 
Chaurasia’s Decoded(dB) 30.97 26.32 33.03 31.06 28.87 x 
[5] Computations(%) 100 100 100 100 100 x 
Predicted(dB) 30.96 26.25 32.90 30.94 28.77 x 
Previous’ Deviation(dB) 0.01 0.07 0.13 0.12 0.10 0.09 
å Computations(%) 68.97 73.05 71.34 51.61 49.95 62.98 
Predicted(dB) 30.97 26.27 32.92 30.94 28.80 x 
Proposed Deviation(dB) 0.00 0.05 0.11 0.12 0.07 0.07V 
Computations(%) 60.96 66.67 63.94 45.12 45.61 56.46 Y 
Chaurasia’s Decoded(dB) 26.56 22.43 28.61 27.17 25.11 x 
[5] Computations(%) 100 100 100 100 100 x 
Predicted(dB) 26.46 22.17 28.42 27.04 24.89 x 
Previous’ — Deviation(dB) 0.10 0.26 0.19 0.13 0.22 0.18 
g Computations(%) 66.89 74.12 67.97 55.76 49.80 62.91 
Predicted(dB) 26.51 22.23 28.48 27.06 25.02 x 
Proposed Deviation(dB) 0.05 0.20 0.13 0.11 0.09 0.09 Y 


Computations(%) 55.18 66.80 59.28 48.73 43.26 54.65 V 


Tab.3 Performance comparison between the previous and proposed methods for Zheng’s method", 


Test images Bird Owl Trees City Cat Average 
Decoded(dB) 30.86 26.52 33.04 31.03 28.94 x 
Zheng’s!°! 
Computations(%) 100 100 100 100 100 x 
Predicted(dB) 30.92 26.48 33.11 31.08 28.89 x 
Previous!“ Deviation(dB) 0.06 0.04 0.07 0.05 0.05 0.05 
4 Computations(%) 79.87 90.77 76.59 55.69 51.92 70.97 
Predicted(dB) 30.92 26.49 33.13 31.08 28.91 x 
Proposed Deviation(dB) 0.06 0.03 0.09 0.05 0.03 0.05= 
Computations(%) 72.54 87.10 67.64 48.05 47.86 64.64 Y 
Decoded(dB) 26.43 22.27 28.45 27.02 24.97 x 
Zheng’s!°! 
Computations(%) 100 100 100 100 100 x 
Predicted(dB) 26.20 22.07 28.31 26.89 24.79 x 
Previous’ — Deviation(dB) 0.23 0.20 0.14 0.13 0.18 0.18 
8 Computations(%) 79.11 89.58 76.88 68.17 53.70 73.49 
Predicted(dB) 26.27 22.16 28.38 26.92 24.92 x 
Proposed Deviation(dB) 0.16 0.11 0.07 0.10 0.05 0.10 Y 


Computations(%) 63.90 83.09 67.37 59.61 46.88 64.17 Y 


Tab.4 Performance comparison between the previous and proposed methods for Gupta’s method”. 
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Test images Bird Owl Trees City Cat Average 


Decoded(dB) 31.00 26.48 33.21 31.28 28.93 x 
Gupta’s!7] ; 
Computations(%) 100 100 100 100 100 x 
Predicted(dB) 30.99 26.35 33.06 31.24 28.84 x 
Previous!**1 Deviation(dB) 0.01 0.13 0.15 0.04 0.09 0.08 
4 Computations(%) 69.48 74.00 72.39 51.71 50.12 63.54 
Predicted(dB) 31.01 26.38 33.11 31.25 28.88 x 
Proposed Deviation(dB) 0.01 0.10 0.10 0.03 0.05 0.06 Y 
Computations(%) 61.60 67.90 64.23 45.26 45.85 56.97 Y 
Decoded(dB) 26.64 22.44 28.64 27.24 25.22 x 
Gupta’sl”] f 
Computations(%) 100 100 100 100 100 x 
Predicted(dB) 26.50 22.19 28.46 27.17 24.99 x 
Previous!**1 Deviation(dB) 0.14 0.25 0.18 0.07 0.23 0.17 
9 Computations(%) 67.19 74.22 67.97 56.15 50.00 63.11 
Predicted(dB) 26.56 22.26 28.52 27.19 25.11 x 
Proposed Deviation(dB) 0.08 0.18 0.12 0.05 0.11 0.11V 


Computations(%) 55.47 66.89 59.28 49.12 43.55 54.86 Y 


5. Discussion 


While the threshold of RE,cxc is set at 0.15, by Eq.(18), we have: 


REcace=0.15 (27) 


Then, we have: 


ACE®=6,13 x ACE" (28) 


While the actual percentage of accumulated collage error (APACE) is defined as 


ACE Coded 


APACE= ACE! + ACRU (29) 


Then, the minimum of APACE, APACE mn , can be defined and calculated by Eqs.(29) 


and (28), respectively, as: 


ACE Coded 
Min ACEC" + ACE Uncoded 


Max 


APACE 


~86% (30) 


In the previous work, the encoding process terminated while APACE, reaches 90%. 


Thus, we know that REcace set at 0.15 needs fewer computations than APACE,,,, set at 
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90%. In other words, while RE yc; reaches 0.15, APACE, is just about 86% and 
smaller than 90%, which implies that in the previous method, we need to encode extra 4% 
range blocks to make APACE,,, reach 90%. Thus, the corresponding computations in the 
proposed method can be are saved. For example, while the size of range blocks is set at 4 x 4, 
the green and red blocks in Fig.3 represent the range blocks encoded in the previous method, 
and the green blocks represent the range blocks encoded in the proposed method. For Bird 
image, we totally have 4,096 range blocks. In the previous method, while APACE min 
reaches 90%, we need to encode 2,826 range blocks which are represented by the red and 
green blocks in Fig.3(a) and comprise 68.99% of total range blocks in Tab.1. In the proposed 
method, while REcacęe reaches 0.15, we only need to encode 2,489 range blocks which are 
represented by the green blocks in Fig.3(a) and comprise 60.77% of total range blocks in 
Tab.1. This implies that the computations of 337 range blocks which are represented by red 
blocks can be saved. While the size of range blocks is set at 8 x 8, similar results as in Fig.3 
were illustrated in Fig.4. For Bird image, we totally have 1,024 range blocks. In the previous 
method, while APACE, reaches 90%, we need to encode 689 range blocks which are 
represented by the red and green blocks in Fig.4(a) and comprise 67.29% of total range 
blocks in Tab.1. In the proposed method, when RE ace reaches 0.15, we only need to 
encode 568 range blocks which are represented by the green blocks in Fig.4(a) and comprise 
55.47% of total range blocks in Tab.1. This implies that the computations of 121 range blocks 
which are represented by red blocks can be saved. Moreover, ACEC in the proposed method 
tries to reflect the proportional relationship among ACE“, ACE;, and ACES), which 


has not been considered in the previous method. Thus, the proposed method can provide 
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more accurate prediction results. In summary, the proposed method can predict the decoded 


image quality more accurately with fewer computations than the previous method. 


(d) City (e) Cat 
Fig.3 While the size of range blocks is 4 x 4, the green and red blocks represent the range blocks encoded 


in the previous method, and the green blocks represent the range blocks encoded in the proposed method. 


(a)Bird (b) Owl (c) Trees 
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(d) City (e) Cat 
Fig.4 While the size of range blocks is 8 x 8, the green and red blocks represent the range blocks encoded 


in the previous method, and the green blocks represent the range blocks encoded in the proposed method. 


6. Conclusion 

An accumulated collage error coefficient based fractal prediction method was proposed in 
this study. Firstly, the definition of ACEC was introduced to describe the proportional 
relationship among the upper bound, lower bound, and actual value of the accumulated 
collage error. Then, the definition and monotonicity of the relative error of ACEC were also 
defined and discussed, respectively. Finally, while the relative error of ACEC reaches a small 
value, the accumulated collage error of all range blocks can be estimated to directly predict 
the decoded image quality. Experimental results verify that the proposed method can provide 
higher prediction accuracy with fewer computations than the previous method. In future work, 
we will attempt to improve the performance of the fractal prediction method in terms of 
prediction accuracy and the number of computations further. 
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